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ABSTRACT 

We present an analysis of archival Chandra and VLA observations of the EO galaxy 
NGC 1399 and the E2 galaxy NGC 4649 in which we investigate cavities in the surround- 
ing X-ray emitting medium caused by the central AGN. We calculate the jet power required 
for the AGN to evacuate these cavities and find values of 8 x 10'*^ erg s~^ and ~ 14 x lO"*^ 
erg s^^ for the lobes of NGC 1399 and ^ 7 x 10*^ erg s^^ and 6 x 10'*^ erg s~^ for those of 
NGC 4649. We also calculate the k/ f values for each cavity, where k is the ratio of the total 
particle energy to that of electrons radiating in the range of 10 MHz to 10 GHz, and / is the 
volume filling factor of the plasma in the cavity. We find that the values of k/ f for the lobes 
of NGC 1399 are - 93 and - 190, and those of the lobes of NGC 4649 are - 15000 and 
^ 12000. We conclude that the assumed spectrum describes the electron distribution in the 
lobes of NGC 1 399 reasonably well, and that there are few entrained particles. For NGC 4649, 
either there are many entrained particles or the model spectrum does not accurately describe 
the population of electrons. 

Key words: galaxies: active — galaxies: nuclei — galaxies: individual: NGC 1399 — galax- 
ies: individual: NGC 4649 — radio continuum: galaxies — X-rays: galaxies 



1 INTRODUCTION 

Within nearby galaxies, groups, and clusters, embedded active 
galactic nuclei (AGN) have been found to interact with the sur- 
rounding hot, X-ray emitting thermal gas, c ausing such distur- 
' bances as s hocks, ripples, and cavities (e.g. Fabian et ajj 1200(1 
' l200ll2006l : lFormanet al.ll2005h . The cavities, which appear as X- 
ray surface brightness depressions, consist of low density relativis- 
tic plasma and are a consequence of the thermal gas being displaced 
by the jets of the AGN. They have been found to be a ssociated with 
the central AGN of many clusters, suc h as Perseus ( Bohringer et al. 
19931 : iFabian et alj|200(]ll200ll2006l). Hydra A (McNamara et al.. 



20od) . M 87 dChurazov et alJl200ll:lFonnan et alJboOSh . and CenT 
taunis jTavlor et al.ll2002l . l20o'^ AGN-induced disturbances are 
also present within the hot int erstellar medium (ISM ) in the haloes 
of normal elliptical galaxies i Diehl & Sta tler 2006). Many early- 
type galaxies, such as M 84 jFinoguenov & Jones,.2001.) . have been 
found to harbor cavities in their X-ray haloes. Our interest in these 
AGN-induced cavities stems from the fact that they can be used 
as calorimeters. The kinetic energy of the jets emanating from the 
central black hole can be estimated by calcula ting the amount of 
energy required to inflate the observed cav ities jSirzan et alj2004l : 
iRaffertv et al.l2006l : lDunn & Fabiaj200^ . Using this estimate, the 
jet power can be calculated once an approximate age of the cavity 
has been deduced. 

For nearby large, elliptical galaxies, accurate measurements of 
the density and temperature profiles of the thermal gas can be ob- 



tained from high resolution Chandra observations. For those sys- 
tems that are sufficiently close, such that radii within an order of 
magnitude of the accretion radius are resolved, these measurements 
can be used to calcu late the Bond i accretion rate, given estimates of 
the black hole mass. lAUen et al.l f2006), in a study of 9 X-ray lumi- 
nous sources that harbor cavities, found that a tight physical corre- 
lation exists between the power available from Bondi accretion of 
the hot gas and the observed jet power. These results are significant, 
especially for studies of accretion and jet formation, as well as the 
formation and ev olution of galax ies. In particular, building upon 
earher studies by jPi Matteo et all ( l2003h and [Taylor et al.i ilOO^ . 
they confirmed that the bulk of the energy produced by the c entral 
AGN in these systems, which are all lFanaroff & RilevI ( 1 19741) type- 
I sources, is released in the relativistic jets. The energy injected 
back into the environment is in principle sufficient to balance ra- 
diative cooling of the hot X-ray gas. T his possibility has also been 
investigated bv lBest et al. Il l2005ll2006i) . who showed that the cool- 
ing of hot gas from the atmosphere of the host galaxy can feasibly 
provide the fuel for low power, radio-loud AGN, and that the heat- 
ing by AGN feedback can balance this cooling, except perhaps in 
the most massive clusters. However, the manner in which the en- 
ergy provided by the AGN is coupled to the thermal energy of the 
hot gas in order to offset cooling, as well as the efficiency of this 
process, are as yet unknown. 

Because of the significance of AGN activity, it is desirable to 
investigate further the population of galaxies which exhibit inter- 
actions between the central AGN and the surrounding medium. In 
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particular, we are interested in those sources which harbor cavities, 
as they offer insight into the physical properties of the AGN itself. 

We analyzed Chandra and VL^B data for NGC 1399 and 
NGC 4649, both of which are large, elliptical galaxies that harbor 
cavities and are located in nearby clusters. Because of their proxim- 
ity, we are able to determine the density of the X-ray emitting gas 
down to radii within a factor of 10 of the accretion radius, making 
these two galaxies ideal for future analysis concerning Bondi ac- 
cretion. At present, we use the cavities in these systems to estimate 
the kinetic energy of the jets and jet power. We also investigate the 
particle content of the cavities by determining fc//, where k is the 
ratio of the total particle energy contained in the cavity to the en- 
ergy accounted for by electrons emitting synchrotron radiation in 
the range of 10 MHz to 10 GHz, and / is the volume filling factor 
of the relativistic plasma in the cavity. The cavities are filled with 
relativistic charged particles and magnetic fields, which together 
determine the synchrotron emission from the lobes. The relative 
energy densities of these two components, and hence k and /, can- 
not be disentangled from the emission alone. We must therefore 
consider the ratio k/ f, which, under the assimiption that the gas 
pressure is equal to the total pressure from the radio plasma, can be 
related to the magnetic field strength. The maximum field strength 
is obtained for the minimum value of k/f — 1, and would indicate 
a pure electron-positron plasma uniformly filling the cavity. 

Our discussion will proceed as follows. We begin in [|2]by de- 
scribing the radio data selection and processing. In ^ we describe 
the X-ray data and preparation, as well as the identification of cav- 
ities. In 21 we present Chandra/WhA comparisons for NGC 1399 
and NGC 4649 and determine the power of the jets in these sys- 
tems. We present the k/f ratios in ij5] followed by conclusions in 

m 

Throughout this discussion, we assume a flat ACDM cosmol- 
ogy with = 0.3 and i^o = 70 km s"^ Mpc"\ 



2 RADIO DATA: REDUCTION AND IMAGING 

The VLA radio data were obtained from the NRAO archive. In 
selecting the data to be analyzed, preference was given to obser- 
vations with time on source J5 5 minutes. Because steep spectrum 
emission is brighter at lower frequencies, preference was also given 
to observations performed at 1.4 GHz. In addition, A-configuration 
or B-configuration observations at 1.4 GHz were desirable in order 
to provide arcsecond resolution for comparison with X-ray images 
from Chandra. Parameters for the VLA observations of NGC 1399 
and NGC 4649 are tabulated in Table[Tl 

The radio data were reduced in the standard manner using 
AIPS. After an initial editing of the data, absolute amplitude and 
phase calibration were performed on each dataset using the scripts 
VLAPROCS and VLARUN. As the flux calibrator for each dataset 
was resolved, a model was used for the calibration. If bad data were 
still present after the initial calibration, those data were flagged and 
the calibration was repeated. 

Each source was imaged using the task IMAGR. Sidelobes 
from outlying sources were removed by using multiple facets while 
imaging. Proper placement of the facets was determined using the 
task SETFC, which was set up to search a 0.5 degree radius for 



sources in the NVS^I catalogue with flux J? 10 mjy. Typically, 
sources outside of that range are either not bright enough or too far 
from the pointing center to have an appreciable effect on the quality 
of the image. As both sources had sufficient signal-to-noise ratios, 
imaging and phase-only self-calibration were then performed iter- 
atively, until the theoretical noise was reached or until the quality 
of the map ceased to benefit from the iterations. Radio flux density 
measurements are shown in Table|2l 



3 X-RAY DATA PREPARATION 

The Chandra X-ray data were used to determine the deprojected 
temperature and density profiles of the X-ray emitting gas in the 
galaxies. These profiles, in combination with the measured sizes of 
the cavities identified from the X-ray and radio emission, allow the 
P dV work done by the cavity on the surrounding gas halo to be 
calculated. The details of the X-ray observations are summarized 
in Table[T] 

Before analysis, the X-ray data were reprocessed and cleaned 
using the CIAO software and calibration files (CIAO v3.3, CALDB 
v3.2). We began the reprocessing by removing the afterglow de- 
tection and re-identifying the hot pixels and cosmic ray afterglows, 
followed by the tool ACIS_PROCESS_EVENTS to remove the pixel 
randomization and to flag potential background events for data 
observed in Very Faint (VF) mode. The Charge-Transfer Ineffi- 
ciency (CTI) was corrected for, followed by standard grade selec- 
tion. Point sources were identified using the WAVDETECT wavelet- 
transform procedure. As these two galaxies were observed with the 
ACIS-S3 chip, background lightcurves to check for flaring were 
taken from the ACIS-Sl chip. For the spectral analysis, back- 
grounds were taken from the CALDB blank-field datasets. They 
had the same reprocessing applied, and were reprojected to the cor- 
rect orientation. 

Spectra were extracted in annular regions centered on the 
peaks of the X-ray emission from the galaxies. Annular regions 
were automatically assigned with constant signal-to-noise, stop- 
ping where the background subtracted surface brightness became 
consistent with zero. The initial signal-to-noise was 100, and this 
was increased or decreased by successive factors of \/2 to obtain a 
number of regions between four and ten. The minimum signal-to- 
noise allowed was 10. 

The 0.5 — 7keV spectra were extracted, binned with a 
mini mum of 20 c ounts per bin, and, using XSPEC (vl2.3.0) 
("e.g. lAmaudlll996h . a PROJCT single temperature MEKAL (e.g. 
iMewe et~ll995h model with a PHABS absorption was used to de- 
project the cluster. The deprojected cluster temperature, abundance 
and normalization profiles allowed the calculation of density and 
pressure profiles. Chi-square was used as the estimator in the fits. 

These profiles give azimuthally averaged values for the cluster 
properties and have been used in the subsequent calculations. The 
substructure present in these galaxies is not expected to bias the 
results on the gas properties significantly ( Donahue et al. 2006). As 
such the use of these azimuthally averaged values is not likely to 
introduce large biases into the subsequent calculations. 



^ The VLA (Very Large Array) is operated by the National Radio Astron- 
omy Observatory (NRAO). 



2 The NVSS (NRAO VLA Sky Survey) is a 1. 4 GHz continuum s urvey 
covering the entire sky north of -40° declination iCondon et all 19981) . 
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Table 1. Positions, redshifts, and radio and X-ray observation details for NGC 1399 and NGC 4649. The frequency (in GHz), configuration, date, and duration 
(in s) of the VLA observations are shown in columns 5 through 8. The ObsID and final exposure time (in ks) of the Chandra observations are shown in columns 
9 and 10. 















VLA 




Chandra 


Source 




Dec^ 


Redshift^ 


Frequency 


Config 


Date 


Duration 


ObsID 


Exposure^ 




(J2000) 


(J2000) 




(GHz) 






(sec) 




(ks) 


NGC 1399 


03h 38m 29.08s 


-35d 27m 02.7s 


0.004753 


1.465 


H3 


1983 Dec 16 


27910 


319 


55.9 


NGC 4649 


12h43m 39.66s 


+1 Id 33m 09.4s 


0.003726 


1.465 


B 


1984 Jan 24 


4410 


785 


36.9 



^Positions and redshifts obtained from the NASA/IPAC Extragalactic Database (NED). 
^The exposure after reprocessing. 
''Hybrid configuration (A and B). 



Table 2. Radio flux density measurements 



Source 


Peak Flux 


Total Flux Density 


RMS 


Beam 


PA 




(Jy beam^^) 


(Jy) 


(/iJy beam-i) 


(arcsec^ ) 


(deg) 


NGC 1399 


1.89x10^2 


4.63x10-1 


72.1 


4.15x2.80 


38.5 


NGC 4649 


1.72x10-2 


2.82x10-2 


29.7 


4.51x3.64 


44.8 



Table 3. Properties 





NGC 1399 


NGC 4649 


Distance (Mpc) 


19.4 


17.1 


Lx (erg s-i)* 


4.27 X 10"! 


1.91 X 10*1 


Pi. 4 (ergs-i Hz-i)t 


2.09 X 1029 


9.87 X 1027 


Scale (kpc arcmin— ^ ) 


5.64 


4.96 


Largest angular size (arcsec)t 


233 


46.5 


Largest physical size (kpc)t 


21.9 


3.85 



* X-ray luminosity, from lO'SuUivan et aljj200ll) . 

tTotal isotropic radio power at 1.4 GHz: Pi. 4 = AttD^Sia- 

t Radio source size. 



4 INTERACTION BETWEEN AGN JETS AND 
THERMAL GAS 

In order to investigate any possible interactions between the AGN 
and the surrounding thermal gas, we superimposed the VLA ra- 
dio contours onto the Chandra X-ray images. These X-ray/radio 
overlays are shown in Figures [T] and [2] Panels (a) and (b) of each 
figure show an adaptively smoothed X-ray image and an unsharp- 
masked X-ray image, respectively, with color scales chosen to high- 
light faint features. The adaptive smoothing was don e with a 2 arc- 
sec k ernel up to a maximum of 20 using ASMOOTH jEbeling et al.l 
l2006h . The unsharp-masked images are difference images made by 
first smoothing the raw image using a large Gaussian kernel (80 
arcsec for NGC 1399 and 20 arcsec for NGC 4649), and then sub- 
tracting that image from one which has been smoothed using a 4 
arcsec Gaussian kernel. The radio contours aid in the determination 
of the approximate size and location of the cavities. We use these 
estimates, in conjunction with density and temperature profiles de- 
rived from the X-ray data, to estimate the jet power required to 
inflate the observed cavities. These calculations are summarized in 
^43] following individual discussions of NGC 1 399 and NGC 4649 
in ij4. 1 1 and ^4.2\ respectively. Some additional derived properties 
of the galaxies are presented in Table[3] 



4.1 NGC 1399 

NGC 1399 is the central EO galaxy of the Fornax cluster. It is host 
to a low powe r radio source, whose lobes are confined within the 
optical galaxy. IPaolillo et all ( l2002h present a ROSAT/VLA com- 
parison which shows interactions between the central AGN and the 
surrounding hot gas. Figure [T] shows our ChandrafVhA compari- 
so n for NGC 1399. The se comparisons confirm the structures found 
m IPaolillo et al] ( |2002|) . as well as the presence of cavities associ- 
ated with the radio lobes. The lobes extend north and south along an 
axis oriented ~10 degrees west of north. They are approximately 
symmetric, each having a length of roughly 1 1 kpc. The central X- 
ray source is elongated in the direction of the lobes, with bulges 
along the east and west edge of the northern and southern lobes, 
respectively (see Figure[TJ)). The southern lobe bends sharply west 
before gradually bending eastward. As can be seen in Figure[TJa), a 
bright rim of X-ray emission follows the western edge of the south- 
ern lobe, marking the edge of the cavity. The northern lobe also 
bends slightly eastward at the end, and a bright rim of emission can 
be seen in Figure [TJa) to follow the western edge of this lobe as 
well. Both lobes end in a bright clump of X-ray emission. 



4.2 NGC 4649 

NGC 4649, also known as M 60, is an E2 galaxy located in a group 
at the eastern edge of the Virgo cluster. Figure |2] shows our Chan- 
drafVhA comparison for NGC 4649. Figure|2jb) shows cavities in 
the X-ray emission that are coincident with the lobes of the weak ra- 
dio source. The lobes extend northeast and southwest, along an axis 
approximately 31 degrees east of north. The northern lobe bends 
westward, while the southern lobe bends eastward, giving the radio 
source an overall S shape. The central X-ray source can be seen in 
Figure|2jb) to have slight bulges along the axis of the radio source. 
These bulges appear to also be coincident with the axis of the 5 
arcsec long central b ar found by [ Ran dall et al. (2004) in an anal- 
ysis of Chandra data. iRandall et anT2004, ,2006 ) also suggest the 
existence of "fingers" in the X-ray emission which extend radially 
from the center of the galaxy. These features are on a much larger 
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Figure 1. X-ray images of NGC 1399 with 1.4 GHz radio contours. The left panel shows an adaptively smoothed X-ray image with a logarithmic color scale, 
and the right panel shows an unsharp-masked X-ray image with a hnear color scale. The unsharp-masking was done by subtracting a raw image smoothed 
using an 80 arcsec Gaussian kernel from one smoothed using a 4 arcsec Gaussian kernel. The dark ring in the unsharp-masked image is due to the central 
source being much brighter than the rest of the image. 




(a) (b) 

Figure 2. X-ray images of NGC 4649 with 1.4 GHz radio contours. The left panel shows an adaptively smoothed X-ray image with a logarithmic color scale, 
and the right panel shows an unsharp-masked X-ray image with a linear color scale. The unsharp-masking was done by subtracting a raw image smoothed 
using a 20 arcsec Gaussian kernel from one smoothed using a 4 arcsec Gaussian kernel. 
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scale than the radio emission, however, and thus are not present 
in the region shown in our images. We note that there are several 
small radio structures at the lowest contour level on the east and 
west sides of the X-ray source, some of which appear to coincide 
with depressions in the X-ray images. These are likely to be arti- 
facts, though future observations may aid in determining whether 
or not they are real features. 

4.3 Jet Power 

In order to estimate the jet power required to evacuate the observed 
cavities in the X-ray ernitting medium, we employ the method de- 
tailed in l Allen et alj|2006h . Using the radio contours as a guide, we 
model the cavities as ellipsoids with semi-axis dimensions ri along 
the direction of the jet, across it, and along the line of sight. 
As the actual depth of the cavities is unknown, we assume sym- 
metry about the axis defined by the direction of the jet, such that 
ta = r„ . Thus the volume of the cavity is given hyV — A-Kr\r% / 3, 
though systematic uncertainties in the cavity dimensions cause this 
value to be uncertain by about a factor of 2. This uncertainty in 
the volume is taken into account in subsequent calculations. The 
thermal pressure P at the center of the cavity is determined from 
temperature and density measurements of the X-ray emitting gas. 
From these values, we can calculate the PAV work done on the 
X-ray gas by the AGN jets in inflating the cavities. The energy E 
required for this process is then found by adding the PAV work 
done on the X-ray gas to the internal energy of the cavity: 

E^PAV + — ^—PV=—^PV. (1) 



72 — 1 72 — 1 

Here 72 is the mean adiabatic index of the relativistic plasma con- 
tained in the cavity and has a value of 4/3. Thus the energy of the 
cavity is given hy E = APV. 

The jet power can be estimated using the kinetic energy E of 
the jet and the age of the cavity. We approximate the age of the cav- 
ity as the sound-speed expansion timescale tea, which is the time 
required for the cavity to expand to its obse rved size at the speed of 
sound in its current env ironment (see e.g. iDunn. Fabian & Tavloj 
l2005l : lBirzanetalJl2004h : 



■ n/cs, 



(2) 



where Cs is the local sound speed. The jet power can then be esti- 
mated as 



' jet, cs 



= E/tr 



(3) 



Because the cavities consist of low density plasma and are 
thus expected to rise buoyantly through the surrounding medium, 
the age of the cavity can also be approximated as the buoyancy 
rise time tbuoy, which is the time required for the cavity to rise 
buoyantly from the center of the galaxy to its current location. The 
buoyancy velocity is «buoy = ^JlgV j SCt>, where S = 7rr^ is 
the cross-sectional area of the cavity, Cd ^ 0.75 is the drag coef- 
ficient, and the gravitational acceleration is given by g = GM{< 
R)/R^. The age, then, is approximately 

^buoy R/Vhuoy, (4) 

where R is the distance from the center of the black hole to the cen- 
ter of the cavity. The mass profiles used in the calculation of t^buoy 
come from the deprojection of the X-ray data assuming hydrostatic 
equilibrium. The jet power can be calculated as 

P'}cX. buoy -£'/^buoy (5) 

The Pjct calculations are summarized in Table]?] 



4.4 A Note on the Timescales 

We have considered two different timescales in the calculations of 
jet power: one derived from the local sound speed and the other de- 
rived from the buoyancy velocity. The sound speed timescale mea- 
sures the expansion of the cavity at its current location, and is an 
appropriate estimate in cases in which the cavities are active, mean- 
ing that they are currently powered by the jets of the AGN. The 
buoyancy timescale, however, measures the rise time of the cavity 
to its current location, and is an appropriate estimate in the case 
of cavities that have detached from the core and are rising buoy- 
antly through the surrounding medium. Such is the case for ghost 
cavities, which are not associated with GHz radio emission. 

In the case of NGC 1399, tbuoy is on average a factor of 2 
greater than tea- The cavities are fairly well detached from the 
core, as can be seen in Figure[T] However, they are associated with 
1.4 GHz emission and thus are presently powered by the AGN jets. 
Because the cavities in NGC 1399 are detached but still active, it is 
difficult to determine which timescale is the most appropriate de- 
scription for the ages of the cavities in this system. In NGC 4649, 
the resolution of the radio image does not allow us to determine if 
the cavities are detached or not, though both timescales are gener- 
ally in agreement with each other within the errors. 



5 PARTICLE ENERGIES AND VOLUME FILLING 
FRACTIONS 

Using the synchrotron emission from the radio lobes, the en- 
ergy contained within the relativistic electrons can be calcu- 
lated. Under the assumption that the radio lobes are in pres- 
sure balance wi th their surro undings, their particle content can 
be investigated. iFabian et alj l [2002f) stu died the lobes o f 3C 8 4 
in the Perseus cluster, and su bsequently iDunn & Fabiarj ( l2004h : 
iDunn. Fabian & Tavloj i2005l) investigated a large number of radio 
sources in galaxy clusters. For a m ore detailed description of the 
method, see lDunn & FabianI ( [2004h . Here we outline the important 
details and any changes to that method. 

To calculate Ec, the energy in the relativistic electrons, we 
assume a continuous synchrotron spectrum with a single spectral 
index a, defined from Si, ot , between v\ = 10 MHz and V2 = 
10 GHz: 



Ee ^^-K X 10' 



a + 0.5 



-B 



'3/2 



(6) 

(e.g. iFabian et alj|2002l) . We adopt spectral index values of a = 
—0.6 and a = —0.7 for the northern and southern lobes of 



NGC 4649, respective lv jStanger & Warwic k 1986), and a = 
-0.92 for NGC 1399 fekers et alJI 19891) . Also taking into account 
the energy within the magnetic field, the total energy in the lobes is 
given by 



o2 



(7) 



where k accounts for any other particles present in the lobe which 
are not accounted for by the simplistic model spectrum, V is the 
volume of the cavity, and / is the volume filling fraction of the 
relativistic plasma. 

The magnetic field present inside the cavity is estimated by 
comparing the synchrotron cooling time of the plasma to the age 
of the lobe. The latter can be estimated from the sound-speed ex- 
pansion timescale or the buoyancy rise time. Under the assumption 
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Table 4. Summary of the calculation of jet power Pj^t (in units of 10*^ erg s~^). "Lobe" refers to the lobe for which the calculation is done (N and S for 
northern and southern lobe, respectively). R is the distance (in kpc) from the center of the black hole to the center of the cavity, P is the thermal pressure (in 
eV cm~'^) of the X-ray gas at radius _R, V is the volume (in cm'') of the cavity, E is the kinetic energy (in units of 10^* erg) of the jets, tea is the sound speed 
expansion timescale (in 10® yr), and tbuoy is the buoyancy rise time (in 10® yr). 



Source 


Lobe 


R 


P 


V 


E 


^CS 


^buoy 


-fjct, cs 


^ct, buoy 






(kpc) 


(eV cm"'') 


(10®3 cm^) 


(1054 erg) 


(10® yr) 


(10® yr) 


(lO^i erg s-i) 


(10*1 ergs-i) 


NGC1399 


N 


8.4 ±0.8 


16.0 ±2.64 


1300 ± 535 


131 ± 59.2 


5.08 ±0.55 


12.6 ± 3.83 


8.24 ±3.70 


3.24 ± 1.49 




S 


8.4 ±0.8 


16.0 ±2.65 


2950 ± 1210 


300 ± 134 


6.98 ±0.76 


10.8 ± 3.30 


13.7 ±6.09 


8.57 ±3.93 


NGC4649 


N 


1.1 ± 0.1 


133 ± 28.3 


50.0 ±20.4 


41.8 ±19.7 


1.91 ±0.20 


1.32 ± 0.36 


6.95 ±3.26 


9.95 ±4.54 




S 


1.3 ±0.1 


134 ± 28.3 


43.7 ± 18.0 


36.6 ±17.5 


1.94 ±0.20 


1.72 ± 0.48 


6.03 ±2.84 


6.63 ±3.11 



that the lobes are in pressure balance with their surroundings, we 
can obtain upper limits on the ratio k/ f: 

,/ \ Sit J a 

If k/f — 1, the lobe is uniformly filled with a purely electron- 
positron plasma with energies corresponding to emission only in 
the range of 10 MHz to 10 GHz. If / ~ 1 and k/f > 1, then 
fc > 1, which implies that there are "extra" particles present. These 
could be thermal protons mixed into the relativistic plasma as the jet 
travels out from the AGN, or they could be electrons which radiate 
outside of the assumed spectrum. 

We calculate the k/f values using timescales derived from 
the sound speed, as well as those derived from the buoyancy veloc- 
ity. These values are presented in Table |5] The equipartition values 
k/ fcq are calculated under the assumption that the plasma within 
the cavity is at equipartition with the pressure from the surrounding 
medium. We use a Monte-Carlo algorithm to calculate the uncer- 
tainties in the upper limits obtained on k/f. The resulting distri- 
bution of limits on A:// is non-gaussian. We find that a log-normal 
distribution is a fair description of the data, and it is from this that 
the values for the uncertainties quoted in Table [5] are taken. The 
values for the uncertainties are large, approximately a factor of 3 of 
the k/f values. 

There is a striking difference between the k/f limits, namely 
that those of NGC 1399 are much lower than those of NGC4649. 
This difference is linked to the radio flux, which is roughly a factor 
of 20 lower in NGC 4649. The lower the value of k/f, the closer 
the assumed electron energy spectrum is to the actual one in the 
lobe. However, as the synchrotron plasma ages, the high energy 
electrons lose their energy more rapidly, and the distribution of 
energies shifts to lower energies and also spreads out. Therefore, 
observing old radio lobes at high frequencies results in low radio 
fluxes, which, if the spectrum is assumed to be flat, underestimates 
the number of particles present. As a result, there would be a mis- 
match between the internal and external pressures. In these calcu- 
lations we force the lobes to be in pressure balance by increasing k 
and/or decreasing /. 

We currently have a lower cut-off to the synchrotron spectrum 
of = 10 MHz. However, in old lobes there would be very low en- 
ergy electrons, and so this lower cut-off might be too high. Under 
the assumption that k/ f — 1 and that we have correctly modeled 
the synchrotron spectrum, we calculate the appropriate i/i . We set 
an absolute lower limit of the cyclotron frequency appropriate to 
the magnetic field determined within the radio lobes. For the north- 
em and southern lobes in NGC 1399, k/ f — 1 can be obtained us- 
ing a minimum frequency of 265 Hz and 49 Hz respectively (the cy- 
clotron frequencies are 54 Hz and 43 Hz). However, in NGC 4649, 



even with the lower cut off at the cyclotron frequency (~ 100 Hz), 
k/f is 2700 and 890 for the northern and southern lobes respec- 
tively. 

Another explanation is that the assumed spectrum is close to 
what exists in the cavity but that there are non-relativistic parti- 
cles present in the lobes, especially in NGC 4649. Thermal protons 
could be entrained as the jet travels through the inner parts of the 
galaxy and halo. These would exert a pressure on the surrounding 
X-ray gas, but would be undetectable from the radio synchrotron 
emission. 

In NGC 1399, the limits on fc// implies that there are few 
entrained particles and that the assumed spectrum does describe 
the distribution of electrons in the lobe reasonably well. However, 
in NGC 4649 either there are many entrained particles or the model 
spectrum does not describe the electron population well. 

There is close agreement between the limits obtained from the 
sound speed and buoyancy timescales. We find that the only lobes 
which can be in equipartition between the leptons and the magnetic 
field are those from NGC 4649 if the buoyancy timescale is used to 
calculate the magnetic field from the synchrotron cooling time. 



6 CONCLUSIONS 

In analyzing Chandra and VLA data for NGC 4649, we have found 
cavities in the surrounding X-ray emitting gas that have not been 
previously reported. We used these cavities and those in NGC 1399 
to estimate the power of the jets emanating from the central black 
hole of these systems. The Pjot values for NGC 1399 using the 
sound speed expansion timescale are ~ 8 X lO''^ erg s~^ for the 
northern lobe and ~ 14 x 10^^ erg s~^ for the southern lobe. The 
values using the buoyancy rise time are on average a factor of 2 
lower for both lobes, but are roughly in agreement with the afore- 
mentioned values within the errors. As the cavities in NGC 1399 
are detached from the core but presently powered by the AGN jets, 
it is difficult to determine which timescale is the relevant one in this 
case. In NGC 4649, the Pj^t values using the sound speed expan- 
sion timescale are ~ 7 x 10^^ erg s~^ for the northern lobe and 
~ 6 X 10*^ erg s^^ for the southern lobe. Calculations using the 

buoyancy rise time yield similar values, within the errors^ 

Comparing these jet powers to those calculated in lAUen et al.l 

1 I2OO6I) , we find that our sources have jet powers that are lower than 
most of those in their sample. Examples are M 87 and NGC 4696, 
both of which are located within dense cluster environments. M 87 
is near the center of the Virgo cluster and has jet powers which 
are roughly an order of magnitude larger than those of our sources. 
The jet powers for the lobes of NGC 4696, which is in the center of 
Centaurus, are roughly 3 to 4 times those of our sources. We find 
that our sources have similar jet powers to those of NGC 4552 in 
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Table 5. k/ f values. S is the radio flux (in Jy) of the specified lobe, k / fcq is the equipailition value, k / /sound is calculated using the sound-speed expansion 
timescale, and fc//buoy is calculated using the buoyancy rise time. 



Galaxy 


Lobe 
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k/ fcq 


k/ /"sound 


k/ /buoy 


NGC 1399 
NGC 4649 


N 
S 
N 
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0.11 ±0.01 
0.13 ±0.01 
0.0038 ± 0.0002 
0.0034 ± 0.0002 
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Virgo, which has values of ~ 7 x 10'*^ erg s"^ and ~ 8 x lO"*^ 
erg s~^, as well as those of NGC 5846, which dominates a group 
of intermediate mass and has jet powers of ~ 4 x 10*^ erg s~^ for 
each lobe. 

We also investigated the particle content of the cavities by cal- 
culating k/ f. We find that the k/ f values for NGC 4649 are much 
larger than those of NGC 139 9 by about 2 orders of magnitude. 
iDunn, Fabian & Tavloil l l2005l) found that there is a large range of 
values for k/f in cavities found in cluster galaxies. In general, k/ f 
is in the range of ~ 1 to ~ 1000 for cavities that are active. The 
values for NGC 1399 fall into this range; however, in the case of 
NGC 4649, the values for k/f using either timescale are more than 
a factor of 10 larger than the upper limit of what has been previ- 
ously observed. 

The minimum value of fc// — 1 corresponds to a lobe which 
is filled with a uniform electron-positron plasma with energies cor- 
responding to emission in the range of 10 MHz to 10 GHz. Large 
values of fc//, such as those found in NGC 4649, can be achieved 
by large values for k or small values for /. Large values for k imply 
that, in order for there to be pressure balance, there must be parti- 
cles present within the lobe which are not accounted for by the as- 
sumed spectrum. These particles can be electrons whose emission 
is outside the range of 10 MHz to 10 GHz or thermal protons that 
either exist in the jets during their formation or are entrained as the 
cavities travel through the surrounding medium. The filling factor 
/ describes the fraction of the volume of the cavity that is occupied 
by radio emitting plasma. The remaining volume fraction can either 
be occupied by thermal gas or plasma whose relativistic electrons 
have aged. Dunn, Fabian & Taylor (2005), in studying ghost cavi- 
ties as well as active cavities in Centaurus and Perseus, found that 
the older cavities have larger k/f values than the active ones and 
have attributed this to the aging of the relativistic electrons. 

The k/f values for NGC 4649 are similar to those of ghost 
cavities in which the relativistic electrons have aged, though the 
values can also suggest that the the lobes contain entrained parti- 
cles. If the electrons in the lobes of NGC 4649 have aged, the model 
spectrum is simply not an accurate description of the electron distri- 
bution. This scenario seems likely, as the radio flux is so small that 
NGC 4649 is expected to soon become a ghost cavity system should 
clear cavities remain. For NGC 1399, the model spectrum appears 
to be an appropriate description of the lobe contents, and there 
are probably not man y entrained particles. As was proposed by 
iDuiin & Fabiad l2004h , higher power radio sources are less likely 
to pick up extra particles from the surrounding medium. Our results 
agree with this description in that NGC 1399 has lower k/f values 
than NGC 4649 and also higher radio power. However, there are 
many physical parameters to consider in investigating these values, 
and no clear correlation has been found between k/f and any one 
parameter. 

In calculating k/ f we. have assumed a simple power-law spec- 
trum. Using high frequency observations to obtain radio flux mea- 



surements of a lobe whose relativistic electrons have aged would 
tend to cause this power-law approximation to underestimate the 
energy of the relativistic electrons in the lobe and therefore overes- 
timate fc//. Our results, then, could be said to suggest that the lobes 
in NGC 4649 are likely to be older than those in NGC 1399. How- 
ever, these results could also suggest that the power-law approxima- 
tion may be poor in both cases. In order to determine whether the 
power-law assumption is reasonable, it would be desirable to ob- 
tain sensitive lower frequency VLA observations of these galaxies, 
such as at 330 MHz, that would allow for accurate measurements 
of the radio flux in the regions corresponding to the cavities. 
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